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MESS 


(Misalignment Estimation Software System) 
for In-Flight Alignment and Calibration 
of Spacecraft Attitude Sensors 

R. desJardins 

Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

Spacecraft with very high precision slewing and attitude maintenance re- 
quirements typically combine a number of high precision celestial reference 
sensors with one or more inertial reference sensors. It then becomes necessary 
to have the capabilities of aligning these various attitude reference sensors 
relative to one another in flight, and of making in-flight calibrations of the slew 
angle scale factors and drift rates of the gyros. The requirement to do this on 
the OAO led to the system of ground and spacecraft procedures and analysis 
known as MESS. Attitude sensor errors arising from preplanned spacecraft 
activity are returned to the OAO support computer data base. Then the MESS 
subroutine system correlates these errors with the spacecraft activity and 
processes these data to arrive at estimates of the sensor alignment and calibra- 
tion parameter values which caused the errors. MESS then automatically com- 
putes correction factors for the computer model of the spacecraft, so that sub- 
sequent spacecraft commands will reflect compensations for the misalignments 
and calibration deviations detected. 
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1. Introduction 

MESS (Misalignment Estimation Software System) is a system of subroutines, 
data sets and spacecraft and ground systems operational procedures for estima- 
ting the inflight values of certain attitude sensor alignment and calibration pa- 
rameters on orbiting spacecraft. Typically MESS is applicable to spacecraft 
with very high precision slewing and attitude maintenance requirements. Such 
spacecraft will normally combine a number of high precision external reference 
sensors (star trackers, sun sensors, possibly earth sensors) with one or more 
internal reference sensors (inertial reference assemblies). To achieve very 
high precision in slewing and pointing, it becomes essential to be able to align 
these various attitude reference sensors relative to one another in flight, and 
also to calibrate in flight the slew angle scale factors and drift rates of the gyros. 

The version of MESS currently implemented, MESS/OAO, was developed 

r ri> / V v ~ r 

for the Orbiting Astronomical Observatory B (OAO-B). MESS/OA(b has been 
used off-line for some alignment tasks on OAO-A2^and a modified version is 
being developed for OAO-C. OAO-B contained five gimballed star trackers 
(GST), a boresighted star tracker (BST), a fine error sensor (FES) utilizing the 
experiment optics, and an extemely stable inertial reference unit (IRU). How- 
ever, the MESS architecture and mathematical and operations analysis are quite 
general, in the sense that a version of MESS can easily be implemented to ac- 
commodate an arbitrary three-axis stabilized spacecraft, and arbitrary attitude 
sensors aboard that spacecraft. 

For example, MESS/OAO is capable of estimating the inflight values of the 
following parameters: 
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• The misalignments and zenith angles of the gimballed star trackers 
(GST). 

• The misalignments of the boresighted star tracker (BST) and the GEP 
fine error sensor (FES). 

• The slew axis misalignments and slew angle scale factors of the inertial 
reference unit (IRU) . 

• The drift rates of the IRU gyros. 

MESS estimates are all made in conceptually the same way, as follows. If 
all the parameters in question had their nominal values, then under equilibrium 
conditions, there would be no errors in the various attitude sensors. That is, 
each GST would track its star with zero errors, the IRU would slew the space- 
craft to its target attitude with zero errors, etc. In point of fact, however, errors 
do occur in these sensors. Some of these errors are random in nature, but 
systematic errors (biases) are also present in general. For most of these sys- 
tematic errors, the error-causing mechanisms have been identified and modeled 
as first-order corrections to various spacecraft characteristic parameters. 

The MESS system reads the various attitude sensor errors from the telem- 
etry returning from the spacecraft, correlates the errors with the spacecraft 
activity which produced them, and then processes these correlated data to ar- 
rive at best corrected estimates of the parameter values which would cause 
such errors. These estimated parameter values subsequently become the new 
nominal values, putatively producing zero (or at least significantly smaller) 
errors. 
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2. MESS Exercises 


Each MESS capability for aligning or calibrating an attitude sensor corre- 
sponds to a well-defined segment of spacecraft activity, called a MESS Exercise. 
The activity specified for each MESS Exercise is designed to produce the type 
and number of attitude sensor errors which will yield mathematically optimal 
estimates of the underlying sensor alignment and calibration parameters. 

The sequence of ground system activity which occurs before and after a 
MESS Exercise is shown schematically in Figure 1. This activity begins with 
an analysis by the MESS Cadre of the alignment requirements for a particular 
sensor. Using the mathematical model of the spacecraft from the data base, a 
MESS Exercise is generated appropriate for the alignment required. The com- 
mands which will effect the Exercise are then fabricated, and uplinked to the 
spacecraft at the appropriate time. 

After the spacecraft performs each segment of the Exercise, the resulting 
attitude sensor errors are downlinked to a storage disk on the MESS computer. 
When the total Exercise has been completed, the MESS Program processes the 
data stored on the disk to arrive at estimates of the alignment and calibration 
parameter values which they imply. When these values have been validated by 
a MESS Analyst, they are entered into the data base to revise the spacecraft 
model. Subsequent commands fabricated against this revised model should 
then result in improved spacecraft slewing and attitude control. 

Four basic Exercises are used in MESS/OAO. These Exercises are per- 
formed in a particular order during the early mission orbits, so that the results 
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Figure 1. Sequence of Events for MESS Exercise 
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obtained by MESS from any given Exercise can be reflected in the commands 
sent to the spacecraft to effect the next Exercise. 

The DRF Exercise for estimating the drift rates of the IRU gyros results 
in a nominally drift-free gyro platform. The GST Exercise for aligning the 
gimballed star trackers relative to one another results in an internally self- 
consistent set of gimballed stellar reference sensors. Then the PAX Exercise 
is performed to align the pointing axis sensors relative to the consistent GST 
set. Finally, when all stellar sensors have been aligned into a consistent set, 
the IRU Exercise is performed to determine the alignment of each IRU slewing 
axis relative to the stellar reference, and to calibrate the slew angle scale 
factors for precise slew angle sensing. 

2.1. The DRF Exercise. IRU gyro drift causes gradual divergence between 
the celestial reference maintained by the stars in fixed space and the inertial 
platform reference defined by the gyro errors. If it were possible for the 
spacecraft to hold with a combination of stellar sensors on a fixed star pattern 
for a long interval of time (several orbits), the gyro drift would appear as a 
gradual buildup in the gyro error registers. Conversely, if the spacecraft were 
to hold on IRU for a similar time interval, the control system would cause the 
spacecraft to follow the inertial platform reference as it drifted and hence the 
stellar sensors would display a gradual buildup of errors. Of these two possible 
ways of observing gyro drift, the first is computationally more attractive, since 
the error readouts from each channel give the gyro error buildup directly. 

Unfortunately, however, the spacecraft is unable to hold to a fixed stellar 
reference throughout the orbit. For when the stellar pattern changes as a 
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guide star becomes occulted or unocculted, the spacecraft attitude will change 
slightly due to misalignments of the trackers. Thus IRU error readings corre- 
sponding to a second pattern cannot be combined with readings corresponding to 
the first pattern. 

The star patterns do tend to repeat, however, once per orbit. When a pattern 
repeats at two separated time points, the spacecraft can be put under stellar 
control at each, and the IRU errors noted. The change in error for each channel , 
divided by the time interval between readings , gives the basic error buildup rate 
for that channel. If such readings are available at several separated time points 
for the same pattern, good statistical (least-squares) estimates of the rates for 
each channel can be obtained. 

Furthermore, stellar patterns other than the initial pattern can be used to 
provide additional independent estimates of the basic error rates. These data 
are processed in precisely the same way as data from the initial pattern. In 
fact, so long as data taken from any one pattern at one time is combined only 
with data from that same pattern taken at other times, the error readings will 
not be contaminated with spurious IRU platform movement due to pattern changes. 
Thus each pattern which repeats periodically can be used to obtain independent 
estimates of the basic gyro drift rates. 

Against the preceding background discussion, the DRF Exercise philosophy 
can now be presented. An interval of about three orbits is set aside during which 
several different stellar patterns will be available and will repeat cyclically 
during successive station passes. The IRU error registers are reset to zero 
at the beginning of the Exercise, and then not reset again during the Exercise. 
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Thus the IRU reference is allowed to drift. Each time the spacecraft is available 
over a station, it is placed under stellar control at whatever pattern is available 
and allowed to settle out. A telemetry transmission containing attitude sensor 
data is then sent to the MESS/OAO computer. The data of interest in this trans- 
mission are the IRU errors, the time point at which the data are taken, and the 
GST configuration prevailing (gimbal angles and tracker status bits). Such trans- 
missions are simply stored on a disk until the Exercise is completed. 

When the last transmission has been received, MESS DRF processing is 
initiated. The telemetry received is merged with the known spacecraft 

attitude and characteristics and passed to the DRF components of MESS. 

These components automatically determine the stellar pattern which 

prevailed during each transmission, and then sort the transmissions into 
groups, each group corresponding to one pattern. When this has been done, the 

least-squares error rate estimates for each channel, for each pattern independ- 
ently, are derived and printed out, and MESS DRF processing halts. The print 
out, containing the several independent estimates of drift rate in each channel, 
is then handed over to S & C (Stabilization and Control) personnel for interpre- 
tation. These personnel determine the numerical drift rate value to adopt for 
each channel and the factors which must be applied to the gyros to compensate 
for these adopted drift rate values. These compensation factors may then be 
commanded into the spacecraft from the ground during a subsequent contact. 

2.2. The GST Exercise. The gimballed star trackers on OAO in an ideal 
sense constitute angle-measuring devices. If they made their measurements 
perfectly, it would be possible to determine, by mathematical analysis of the 
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gimbal angles, the angle between two given stars being tracked by a given pair 
of star trackers. Due to unknown misalignments of the trackers, however, the 
commanded gimbal angles drive the trackers to a configuration which does not 
exactly match the true angle between their guide stars, and hence when allowed 
to track these stars, the star trackers take on errors caused by the misalignments. 
The mathematical analyst can derive the relationship between the misalignments 
and the corresponding errors. Different stellar configurations will provide 
different equations relating the misalignments to the errors. Given a sufficient 
number of sufficiently different such equations, and given the errors, it is pos- 
sible to solve for the misalignments which produced them. The purpose of the 
GST Exercise is to provide a number of different stellar configurations so that 
the GST misalignments can be determined. 

The GST Exercise is performed as follows. A time interval of at least one 
orbit is set aside during which the spacecraft will make passes over a number 
of different stations, so that different stellar configurations will be available. 

At the beginning of each pass, the spacecraft is put under stellar control at any 
convenient stellar pattern, the IRU is reset, and then the spacecraft is placed 
under IRU control for the remainder of the pass. In general, several of the 
trackers will each now have several guide stars available within their gimbal 
capability. All such trackers are then gimballed together, each to such an 
available star. When the trackers have all acquired their stars and settled out, 
a burst of attitude telemetry data are taken and sent to the MESS computer. 

The data of interest in this transmission are the GST configuration and errors. 
Since the spacecraft is being held on IRU, these errors can be averaged together 
to remove the measurement jitter. The trackers are then all gimballed to 
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different stars and another transmission sent to the MESS computer. Several 
different such configurations can be viewed during each pass, for a number of 
passes at different stations, and a transmission sent to the MESS computer for 
each. Such transmissions are simply stored in a data set until the Exercise is 
completed. 

When the last transmission has been received, MESS GST processing is 
initiated. The telemetry received is merged with the spacecraft attitude and 
characteristics and passed to the GST components of MESS. These components 
compute, for each pair of stars in each configuration, the angle between the stars 
as measured by the misaligned trackers. Coefficients are then computed for 
the linear equation which relates the unknown misalignments of the two trackers 
in each possible pair to the deviation of the measured angle from the true angle 
between that pair of guide stars. Such equations are collected into a linear 
system of equations in the unknown misalignments, and this system is solved in 
the least-squares sense. The resulting estimated misalignments are relative 
misalignments, in the sense that they are nominally consistent among the tracker 
themselves, but still contain an unknown misalignment of the GST set with re- 
spect to the spacecraft control coordinate system. 

At the conclusion of MESS GST processing, the GST misalignments (possibly 
including corrected zenith angles) are entered into the SCPS data base. Subse- 
quent spacecraft commands generated by SCPS will then reflect the improved 
GST alignment parameters. 

2.3. The PAX Exercise. The next task in the OAO attitude sensor align- 
ment sequence is to determine the alignment of the pointing axis sensors - BST 
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and FES - with respect to the self-consistent GST set. This is a fairly simple 
task. If the spacecraft is held on the GST set in all three axes at an attitude 
which nominally points the PAX sensors at a suitable target star, each PAX 
sensor will independently take on errors corresponding directly to its misalign- 
ment. Conversely, if the spacecraft is held on one of the PAX sensors - say 
the FES - in pitch and yaw, and GST in roll, then the GST set will take on gimbal 
errors due to the offset position of the pointing axis caused by the FES misalign- 
ment. In addition, the BST will take on errors corresponding directly to its 
misalignment relative to the FES. 

The PAX Exercise is performed simply by holding on one PAX sensor in 
pitch and yaw and GST or IRU in roll. Then a single attitude data transmission 
is sent to the MESS computer. The data of interest in this transmission are the 
PAX sensor errors and the GST configuration and errors. 

When this transmission has been received, MESS PAX processing is initi- 
ated. The telemetry data are merged with the spacecraft attitude and character- 
istics and passed to the PAX components of MESS. These components compute 
the offset of the GST set from the PAX sensor which is controlling pitch and 
yaw. The pitch and yaw offsets represent directly the misalignments of that 
PAX sensor. The errors of the other PAX sensor give its misalignments directly. 
These misalignments are all then reference transformed to a consistent refer- 
ence. The reference considered nominal by MESS is defined by the FES in 
pitch and yaw, and a least-squares fit to the zenith positions of the side-looking 
GST set in roll. The resulting set of alignment parameters represents a con- 
sistent alignment of all the stellar sensors, and provides stellar referenced 
coordinates within the spacecraft against which the IRU can be aligned. 
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At the conclusion of MESS PAX processing, the alignment parameter values 
determined in the PAX Exercise are entered into the data base, to reflect im- 
proved alignment of the pointing axis sensors in subsequent commands. 

2.4. The IRU Exercise . The last task in the OAO alignment sequence is 
the determination of the IRU slew axes alignments with respect to the stellar 
reference coordinates, and the estimation of the IRU slew angle scale factors. 
The IRU Exercise is the most time-consuming of the MESS/OAO Exercises, 
and requires the most advance planning. The reason for this state of affairs is 
that the IRU slew characteristics are invisible to the spacecraft while the space- 
craft is holding at a fixed attitude. Only when the spacecraft is commanded to 
perform a slew about a given nominal IRU axis does it become apparent that the 
spacecraft has actually slewed about an axis which is slightly misaligned from 
nominal. And only when the spacecraft is commanded to slew through a given 
angle does it become apparent that the spacecraft actually has slewed through 
a slightly different angle, either too large or too small. Thus it takes a slew 
through a significant angle to generate one data point for the IRU Exercise. The 
mathematical minimum number of properly defined slews to determine all the 
IRU parameters is six, but this allows for no statistical redundancy. Twelve 
slews of proper type and size is considered a minimum number for statistical 
purposes, and eighteen slews will give significantly better results, and include 
a margin for operational failures. Since in general only one slew can be per- 
formed per contact, it is necessary to set aside an interval of perhaps ten orbits 
for precise IRU alignment. 

Each IRU axis is aligned separately, as follows. To be definite, consider 
the alignment of the IRU pitch axis. Consider two target stars visible to the 
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FES, separated by an angle of (say) exactly 20°. There exists a spacecraft 
attitude at which the FES is pointing to the first target star, and from which a 
pure pitch slew about the nominal IRU pitch axis of precisely +20° will nominally 
take the FES directly to the second target star. Now if the spacecraft is initially 
positioned at the attitude which will do this, and then commanded to perform an 
IRU pitch of exactly +20°, two things will take place. First of all, the spacecraft 
will slew off to one side or the other of the target, due to any deviation from 
nominal in the alignment of the IRU pitch axis. Secondly, the spacecraft will 
pitch too long or too short, due to any deviation from nominal in the positive 
pitch scale factor. The deviations of the actual spacecraft terminal attitude 
after the slew from that predicted on the basis of nominal IRU pitch slew param- 
eters represents error data which can be used to solve for the parameter vari- 
ations which produced the errors. 

With the above preliminary explanation, and continuing with the +20° pitch 
slew example, the IRU Exercise philosophy can now be presented. In the back 
orbit prior to some contact at which IRU Exercise data is to be taken, assume 
the spacecraft is positioned at a stable attitude, at optimum roll (roll angle 
which gives maximum solar paddle power for specified pointing axis orientation) 
and pointing to the first target star. A roll slew is commanded out of memory 
which will bring the spacecraft to the exact attitude required for the Exercise, 
as explained above. This roll slew is so timed that settling is completed just 
as the contact begins, in order to minimize the time spent away from optimum 
roll. As soon into the contact as the spacecraft is commandable, it is placed 
under stellar control so that it is known to be at the exact planned attitude. 

Then the IRU is reset and the spacecraft is placed under IRU control. All this 
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activity is performed quickly, so as to be completed before the +20° pitch slew 
is commanded out of memory. This command is timed to take place some two 
minutes into the contact. The spacecraft then proceeds to slew, requiring per- 
haps four minutes to complete slew and settling. The spacecraft will terminate 
the slew at an actual attitude which is slightly different from the planned terminal 
attitude, as explained above. The IRU error registers will now all read zero, 
since the IRU will have satisfied its slew command. However the stellar sensors 
will all be reporting some errors, corresponding to the deviation of the actual 
spacecraft attitude from that planned. Next the spacecraft is placed under stel- 
lar control, and allowed to settle out. This places the spacecraft at exactly the 
planned attitude. Thus the IRU will now have taken on some errors, corresponding 
directly to the attitude deviation in roll, pitch and yaw. At this time, an attitude 
data transmission is sent to the MESS computer. 

At each contact of sufficient duration, one such slew can be performed. In 
general, the slews will have been chosen to represent all possible cominations 
of roll, pitch and yaw slews, positive and negative slews, and small, medium 
and large slews. The transmissions received from the various contacts are 
simply stored in a data set until the Execise is completed. 

When the last transmission has been received, MESS IRU processing is 
initiated. The telemetry received is merged with the corresponding spacecraft 
attitudes and characteristics and passed to the IRU components of MESS. These 
components determine the size of each slew commanded as well as the slew 
axis and direction, then compute the coefficients of linear equations relating 
the unknown misalignment and scale factor parameters to the errors observed. 
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All such equations for all the slews represented in the data are collected into a 
linear system of equations, and solved in the least-squares sense for the mis- 
alignments and scale factors. The resulting estimated misalignments and scale 
factors come out already correctly referenced to the spacecraft stellar coordi- 
nate reference. 

At the conclusion of MESS IRU processing, the alignment parameters and 
scale factors determined in the IRU Exercise are entered into the SCPS data 
base, so that subsequent slew commands generated by SCPS will bring the space- 
craft more precisely to the desired attitudes. 

3. MESS System Architecture 

MESS/OAO is a lOOK-instruction component of the 2000K- instruction Sup- 
port Computer Program System (SCPS). SCPS is the software which makes 
flying the OAO possible. The OAO must be continually supplied with fresh 
commands, 24 hours a day, as much just to keep it stable and under control as 
to enable it to collect scientifically useful data. To do this, SCPS takes inputs 
from the mission and scientific staffs, and fabricates contact messages con- 
taining command memory loads which will keep the OAO functioning properly 
until the follow-on memory load is received. This process has been repeated, 
many times a day, for over 2 1/2 years now to keep OAO-A2 alive and well. 

SCPS resides on a dedicated IBM 360/65 system with 1.5M bytes of core. 
MESS/OAO, an off-line component of SCPS, has an overlay tree structure, 
schematically represented in Figure 2, which allows it to fit within a maximum 
of 450K bytes of core. The root segment of MESS primarily consists of the 
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Driver, which invokes paths specified in the control record input stream. The 
overlay segments correspond to the functional areas of Preprocessing, Analysis, 
Solution and Postprocessing. Any MESS run may be interrupted or a previous 
run resumed at any point. Disk common preserves continuity between runs. 

Sequencing through the MESS subroutines is controlled by records placed 
in an input control stream. This stream is normally written on a 2260 display 
device, with card backup. The control stream is also used to make changes to 
MESS System Common between subroutine executions so as to provide for sev- 
eral consecutive runs against the same data, using different run parameters. A 
summary of all significant MESS output is provided at the 2260 display station, 
so that the analyst may have immediate knowledge of results and can direct 
subsequent runs accordingly. An entire MESS alignment episode can be moni- 
tored, controlled and validated by the MESS Analyst from the display station. 

3.1. Preprocessing. There are two MESS preprocessors - PMESS1, the 
telemetry preprocessor, and PMESS2, the mathematical preprocessor. 

3.1.1. PMESS1 , the telemetry preprocessor, selects and merges specified 
data from the data sets containing the actual returned telemetry data, and the 
spacecraft model and activity predicts data. The telemetry data are first edited, 
averaged and weighted, according to options specified in the control stream, and 
converted to engineering units. The spacecraft model and activity predicts data 
are then searched for those data pertaining to the time contained in the telemetry 
data. When a match is obtained, the pertinent spacecraft characteristics and 
predicted attitude are then stripped out of the predicts data and merged with the 
telemetry data. The end result is the temporary storage data set PMESS, which 
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contains a header record of spacecraft characteristics, and additional records 
of merged and edited attitude and telemetry data. 

3.1.2. PMESS2 , the mathematical preprocessor, then reads the PMESS 
data set and pulls out those data pertinent to the type of Exercise - PAX, GST, 

IRU or DRF - which was performed. These data are then transformed to the 
model definitions, and formatted into arrays suitable for direct input to the 
analysis subroutine to be used. 

3.2. Analysis. Next, the pertinent analysis subroutine is called. There is 
a different one of these for each different Exercise. These subroutines contain 
the mathematical analysis of the problem; they all function by analyzing the data 
and computing the elements of arrays which can subsequently be processed by 

a strictly mathematical algorithm to yield the solution estimates directly. 

3.2.1. PAKPRF is used to estimate the IRU gyro drift rates. These rates 
are derived with respect to specific GST patterns. That is, the IRU errors are 
read once or more per orbit while the spacecraft is holding on a particular GST 
pattern. As the IRU platform gradually drifts away from the stellar reference, 
the IRU errors build up, and the rate of buildup is estimated by MESS. In practice, 
various patterns are used as they become available around the orbit, and 
PAKDRF automatically sorts the data by pattern, so as to derive independent 
estimates of rate for each pattern (up to a maximum of five). The DRF Exercise 
is normally performed prior to any of the other MESS Exercises, so as to pro- 
vide a nondrifting spacecraft as required for portions of the subsequent Exercises. 
The output from PAKDRF takes the form of a set of linear equations which can 
subsequently be solved to give the gyro drift rate estimates directly. 
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3.2.2. PAKGST is used to align the gimballed star trackers. These sensors 


are aligned only relatively - that is, they are aligned to one another so as to be 
self-consistent. This Exercise will normally be performed prior to the PAX 
and IRU Exercises, so that the misalignments determined during these latter 
Exercises will be referenced to a self-consistent GST set. The output from 
PAKGST takes the form of a set of linear equations which can subsequently be 
solved to give the raw GST misalignments directly. 

3.2.3. PAKPAX is used to align the pointing axis sensors - the boresighted 
star tracker (BST) and the GEP fine error sensor (FES). These sensors are 
aligned independently of each other to a reference system defined by the GST. 

This Exercise is normally performed prior to the IRU Exercise so that the 
slew axis misalignments determined during that Exercise will be referenced 

to an internally self-consistent set of stellar sensors. PAKPAX contains its 
own solution algorithm, and outputs the raw BST and FES misalignment estimates 
directly. 

3.2.4. PAKIRU is used to estimate the slew axis misalignments and slew 
angle scale factors. That is, if slews in a certain direction are terminating 
consistently off to one side or too long or short, PAKIRU will estimate the 
values of the parameters which describe such deviations. The IRU is aligned 
to a reference system defined by the best available stellar sensors (ideally, the 
FES in pitch and yaw, and GST in roll). The output from PAKIRU takes the form 
of a set of linear equations which can subsequently be solved to give the estimates 
of the raw slew axis misalignments and slew angle scale factor deviations 
directly. 
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3.3. Solution. One solution subroutine is used to solve the linear systems 
constructed by the subroutines PAKGST, PAKIRU and PAKDRF. (PAKPAX 
contains its own solution algorithm.) 

3.3.1. SOLVEQ is a general least-squares algorithm for solving systems of 
up to 300 equations which are over-determined (more linearly independent 
equations than unknowns) in up to 30 unknowns. The algorithm takes as input 
a coefficient matrix C(300, 30) and a right-hand-side vector E(300), as well 
as a row mask MASKR(300) and a column mask MASKC(30) . The algorithm 
solves the equation CX = E where X(30) is the array of solution values, accord- 
ing to a least-squares criterion, restricted by the masks as follows. 

Stated informally, SOLVEQ takes up to 300 equations linear in up to 30 
unknowns, assumes specified values of some selected unknowns, and then solves 
(in the least-squares sense) selected equations in selected remaining unknowns. 
More precisely, the algorithm selects the rows i of C specified by the values 
MASKR(i) = 1 (MASKR will usually be set so that these are the nonzero rows 
of C). Then the algorithm sets solution values for any parameters X(j) specified 
by MASKC(j) = 2. The algorithm then solves in the least-squares sense the 
equations selected by MASKR(i) = 1 for the unknowns X(j) selected by MASKC(j) = 1, 
assuming the values X(j) selected by MASKC(j) = 2. 

The output from SOLVEQ consists of the solution array X. In the event the 
masked system is rank-deficient (fewer linearly independent equations than 
unknowns), SOLVEQ outputs diagnostic information to identify the source of the 
deficiency. 
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3.4. Postprocessing. There are two MESS postprocessors - REFXFM, 


which reference transforms the various raw misalignment estimates to a con- 
sistent reference, and WRBASE, which writes out the revised parameter values 
for the SCPS data base. 

3.4.1. REFXFM contains models of all the various coordinate systems to 
which the misalignments could conceivably be referenced. References for roll, 
pitch, yaw and the IRU are selected by input parameters IROLL, IP ITCH, IYAW 
and IIRU. Misalignments in roll, pitch and yaw can be referenced to any appro- 
priate individual attitude sensor, or to an appropriate least-squares combination 
of GST. The IRU misalignments can be further separated into sets of orthogonal 
and non-orthogonal misalignments, or combined. The output from REFXFM is 

a consistent set of alignment parameter estimates for all the sensors involved 
in the preceding Exercises. It is in REFXFM that the results of Exercises for 
separate sensors are combined into a consistent set of estimates. 

3.4.2. WRBASE takes the consistent set of misalignment and calibration 
parameter estimates from REFXFM, and combines it with the nominal SCPS 
values (current data base values) to produce a revised set of values for the 
SCPS data base. . The output from WRBASE is a table of values in a form suitable 
for direct entry to members IDL005 and IDL009 of the data base, when validated. 

3.5. MESS Common. All MESS subroutines are (with trivial exceptions) 
argument-free, and pass their inputs and outputs through MESS Systems Com- 
mon. All variables in MESS Common have been assigned unique names, and all 
these names have been, listed in a NAMELIST statement in the MESS Driver. 
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Thus all MESS variables are accessible to the external world via NAMELIST 
read and write statements. 

This capability is used primarily to change the values of selected variables 
between calls to MESS subroutines. Thus it is possible in effect to alter the 
MESS run dynamically during execution, provided the changes desired can be 
predefined. For example, a given set of MESS Exercise data can be processed 
under a series of different options, all within the same execution cycle. The 
first run is made, then the Common variables are changed so as to set up the 
second run, then the second run is made, and so on. The changes to MESS 
Common, as well as the sequence of subroutines to be called, are defined by an 
input stream of control records. 

4. Mathematical Analysis 

A brief sketch of the core analysis is given herein. For more details, 
consult the References. The comprehensive MESS Manual, to be published 
shortly, will be a compendium of all published material on MESS. 

4.1. The Spacecraft Model. For purposes of this discussion, a spacecraft 
in orbit stabilized in three axes will be represented by a right-handed ortho- 
gonal control coordinate system c being maintained at attitude A relative to the 
usual geocentric inertial coordinate system o. That is, direction cosines v of 
any direction fixed in inertial space o are transformed by an attitude matrix 
A into direction cosines Av expressed in control coordinates c. From control 
coordinates, the vector Av must undergo the transformation R g to be expressed 
in the (known) local coordinate system s of attitude sensor s. From this local 
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coordinate system, the vector R s Av must be transformed by the modeled 
misalignments I + dR__ , where I is the identity, to be expressed in the (unknown) 
misaligned coordinate system s'. Figure 3 shows the flow diagram of these 

coordinate systems for a typical sensor s. Additional sensors t, u may 

also be referenced to the control coordinate system, as shown in Figure 3. 



Figure 3. Attitude Sensor Coordinate Transforms 


What is included in R s depends on one’s point of view. For the present 
discussion, R s may be considered the nominal transformation (as determined 
by spacecraft design) , while I + dR s represents the total misalignment. In the 
operational environment, some values for the misalignment parameters will 
already have been determined by preflight measurement or previous inflight 
estimation. In these cases, R s may be taken to represent the misaligned co- 
ordinate system specified by those values of the parameters, while I + dR g 
represents residual misalignments which remain to be calibrated out. 

The exact form of R s depends on the location of sensor S relative to the 
control coordinates c, and will be unique to the particular spacecraft. Trans- 
formations used in MESS/OAO are given in appropriate later sections. 
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The attitude A of the spacecraft is defined by a yaw-pitch-roll Euler angle 
sequence. The control coordinate system axes c are obtained from the inertial 
coordinate system axes o by rotating the latter first by a yaw (positive sense 
about positive z-axis) through right ascension a, next by a negative pitch (nega- 
tive sense about positive y-axis) through declination 8, finally by a roll (positive 
sense about positive x-axis) through roll angle In the matrix formulation, 


A = 


1 

0 

0 \ 

/ cs 

0 

SS 

0 

C/3 

S/3 

| ( 0 

1 

0 

0 

-S/3 

C/3 ) 

\ -SS 

0 

cs 


( Ca So- 

-So. Ca 

0 0 


0 

0 
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/ 


C CS Ca 

-C/3 Sa - S/3 SS Ca 

Sj8 Sa - C/3 SS Ca 


CS Sa 

C/3 Ca - S/3 SS Sa 

-S/3 Ca - C/3 SS Sa 



(C = cos, S = sin), so that x c = Ax o . Conversely, given the attitude matrix A, 
one can determine the right ascension a, declination S and roll angle J3 which 
it represents by the equations 


a = AT AN 2 (k 12 , A n ) 


S = AT AN 



/3 = AT AN 2 (A 23 , A 33 ) 
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4.2. Analysis for the DRF Exercise. The DRF Exercise involves only 
classic least-squares analysis. For each given pattern, suppose n data points 
have been collected. These data have the form (t. , e. ), i = 1, 2, . . ., n, where 
t. is the time the i th data point was collected, and e. is the corresponding array 
of IRU gyro errors in roll, pitch and yaw. The errors are modeled as growing 
linearly with time. Hence, it is assumed that there exist unknown constant arrays 
a, b, such that 

e. = a t ; + b, i = 1, 2, • • • , n 

Expressing these as 3n scalar equations in 6 unknowns, we can write 

e ij = a j t i + bj> i = 1, 2, • • • , n, j = 1, 2, 3 

Each of these equations can be considered as expressing the error e i . as a 
linear function of the unknowns a . , b. : 

e. . = t . a. +1 • b., i = 1, 2, ♦ ♦ • , n, j = 1, 2, 3 

Thus the 3n equations can be expressed in the following matrix form CX = E: 
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In general, the data (t £ , ej are inconsistent, so this equation is overde- 
termined in the parameters a, b. In MESS, the least-squares solution is chosen, 
i. e., the solution vector X is found which minimizes ||CX — E | j. The solution 
parameters a represent the useful MESS output; the solution values of a 1 , a 2 , a 3 
are the gyro error rate estimates in roll, pitch and yaw, respectively. 

4.3. Analysis for the GST Exercise. The six gimballed startrackers on 
OAO-A2 are oriented one at each end of the three coordinate axes, as shown 
in Figure 4. Local tracker coordinate systems x k , y k , z k , k = 1, 2, . . ., 6, are 
defined such that the zenith position of tracker k lies along the +x k -axis, the 
inner gimbal axis at zenith coincides with the y k -axis, and the outer gimbal 
axis coincides with the z k -axis. Within the local tracker k coordinate system, 
the outer and inner gimbal phasing is defined according to the right ascension- 
declination convention: outer gimbal motion about the z k -axis is positive counter 
clockwise, inner gimbal motion about the inner gimbal axis is positive clockwise. 
These relationships are shown in Figure 5. 


Hence a star in the field of view of tracker k with gimbal angles a k , /x k 
(outer, inner, resp.) has local coordinates 
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ZENITH 



OGX = OUTER GIMBAL AXIS 

o- = OUTER GIMBAL ANGLE 

IGX = INNER GIMBAL AXIS 

P = INNER GIMBAL ANGLE 

N = NORMAL TO THE GIMBAL 

PLANE ( the plane determined 
by the two gimbal axes) 


Figure 5. Gimbal Angles 


This transformation is represented by the following coordinate-transformation 
diagram: 



As each of the six local coordinate systems is nominally aligned parallel to 
some control set of axes, the nominal transformation R k from the control axis 
c to the system k has a matrix composed only of 0, del: 



Each such transformation is represented by the following diagram: 



The misalignments modeled were, first, rotational misalignments d^ k , do , 
d0 k (taken positive in the conventional right-hand sense) about each of the 
tracker k coordinate axes x k , y , z , resp. These represent an arbitrary mis- 
alignment of the tracker k gimbal platform relative to the spacecraft structure 
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as a whole. Making the usual small-angle (first-order) approximations 

ca ~ 1, sa £ a ( j a j « 1) 


the misalignments can be represented by a small-angle rotation matrix 


i +<«.„ 




1 d./v -de k 

1 <% 
tlP k -<% 1 


\ 

/ 


Second, there were also modeled shifts of the null position in the inner and 
outer gimbals. A null shift in the outer gimbal cannot be distinguished from a 
misalignment about the tracker z k -axis, since the outer gimbal axis is always 
parallel to the z k -axis. Hence no separate parameter is necessary to represent 
this shift. 


However, a null shift d/3 k in the inner gimbal can be separated from a mis- 
alignment about the tracker y k -axis by taking a large outer gimbal angle, since 
this separates the inner gimbal axis from the y k -axis. (The inner gimbal axis 
is parallel to the y k -axis only when the outer gimbal axis is zero.) 


Because d/3 k masquerades as d(9 k for zero outer gimbal angle, the sense of 
d/3 k has been taken as that of d<9 k , viz., positive in the usual right-hand sense. 
This is opposite to the sense of the inner gimbal angle itself , which is that of 
declination (declination is negative in the usual right-hand sense). 

Thus we have defined the following coordinate transformation models for 
each tracker. In the nominal case, the outer and inner gimbal angles cr k , ^ k , 
resp. , are computed based on the known attitude of the spacecraft and the known 
misalignments df^: 
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These gimbal angles would nominally point the startracker k line-of-sight directly 
at the target star. Due to unknown misalignments dO k , d/3 k , however, and to the 
fact that the true spacecraft attitude could be somewhat removed from nominal, 
the tracker will not in general find the target star at the commanded angles a k , 
fj- k , but rather at measured angles cr' k , fJ-' k slightly different from cr , /i fc . Hence 
including misalignments, we have the following situation: 


(7>3j^Q 1 * d -h » Q 1 * a_ , d fU@ 


+ d® 


K- d 4 


(The negative sign before d/3 k is due to the fact that its sense is opposite to that 
of .) Thus for a given star being tracked by a given tracker, we have the 
following: 



where R' = (I + dfi k ) R k . 

Now suppose the spacecraft is holding at some fixed attitude, and consider 
the tracking startrackers by pairs. On the one hand, the angle between two stars 
individually being tracked by two corresponding startrackers is known precisely 
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from star catalogs. On the other hand, the angle between the tracking startrackers 
as computed from the gimbal measurements will differ from the true angle, and 
this discrepancy is assumed to be due to misalignments of the trackers involved. 
Consider the following diagram: 



The computed dot product between the trackers i and j is given by 

b' = s' c • s? c 

j i 

where s . c are the measured (') coordinates of star i (sj in the control coordi- 
nate system ( c ). From the diagram: 



Similarly for tracker j. 


On the other hand, the true dot product between the stars is known from the 
star catalog, and may be expressed as a function of the (unknown) misalignments, 
and the (known) measured angles: 






Taking the dot product in the i' coordinate system: 




The difference A b = b - b* is a function of the eight unknown misalignments d<£ i# 
diS'j, d <//., d/Sj, d <£,, d0, , d <//. , d/3,, and hence may be expressed to first order in 
differential form: 



/ 3Ab \ 
jo 



d0. 

d0. 

1 


where the differential coefficients are all evaluated in the nominal state, i. e., 
assuming zero misalignments. For any given data reading, the discrepancy Ab 
= b - b' is computed from the known coordinates of the stars and the data (the 
measured angles cr [ ,//. , a'. , jjJ. ). The differential coefficients are also computed 
from the data, as follows: 

Ab =b -b' = sf • sf - s, c • s! c , 


BAb 
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(The coordinate system i" is that of the misaligned startracker i gimbal plat- 
form. Cf. Figure 6). Hence 



Evaluating s j' in the nominal state, this becomes 

/ 0 \ 

- s/4 

\ so 1 / 

Hence 
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In a similar way, 
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At zero misalignments, 


= -i' 

o J 


9 ! \ 

ca i s ^i\ 

S°"i s fj[ 

\ - c < ) 


The coefficients 


'3Ab \ / BAb \ /BAb \ /BAb 

^ o' l 3d9 i Jo' ( 3d *j Jo’ 


can be computed using the above formulas, due to the symmetry of i and j, simply 
by interchanging i and j. Hence the following equation in eight unknowns has 
been generated: 

d<£. 
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dt p. 
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where C . . , C, . are 1x4 matrices: 
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Altogether in the system of the six gimballed startrackers, there are 24 unknowns, 
and the second above equation may be regarded as one equation in 24 unknowns. 
Each pair of tracking startrackers generates one such equation for each data 
reading. A set of three tracking startrackers taken by pairs generates three 
such equations, and in general a set of n tracking startrackers generates (^j 
such equations (some redundant). 

Data readings are collected representing many different values of a, /a for 
all trackers, and the equations described above are generated. In this way a 
large number of equations in the 24 unknown misalignments are generated. This 
system of equations is then solved in the least-squares sense. The solutions 
represent the least-squares estimates sought. 


4.4. Analysis for the PAX Exercise. Consider a specified attitude A 0 at 
which the nominal spacecraft control coordinate system c has its +x c -axis directed 
to a specified FES target. At this attitude, outer and inner gimbal angles , 

i = 1, 2, . . . n, are sent to the spacecraft to point the n tracking GST to their 
respective guide stars. In the spacecraft control coordinate system, then, the 
n guide stars are located by the vectors 



where C = cosine, S = sine, do. is the misalignment matrix 



\ 

/ 




and R. is the matrix relating the nominal tracker i local coordinate system to 
the control coordinate system c. 


When a pointing axis sensor PAX is allowed to assume control of pitch and 
yaw, its misalignments d<9 in pitch and dp in yaw cause the spacecraft to move 
in the directions and magnitudes of the negatives of the above misalignments. 
Thus the spacecraft will move in the positive sense about the +y c -axis through 
the angle -d£, and in the positive sense about the +z c -axis through the angle -dp. 
The spacecraft will thus take up a new attitude A x . At this attitude, the n track- 
ing GST will take on outer and inner gimbal errors Aa-.,A/x. , i = 1 , 2, . . n. 
Thus in the new position of the control coordinate system, the n guide stars are 
located by the vectors 


W 


. =RJ (I 


C (c. + Ao-.) C(/z. +A / u.)\ 
S(a. + Act.) C (/x. + A^) 

SO. +A/j..) J 


, i = 1, 2, • • • , n 


The rotation which the spacecraft undergoes in moving from attitude A 0 to 
attitude A x can be determined by finding the rotation matrix R which will move 
the vectors {v jj- into the vectors {wj in the sense of making R Vj = w ; , 
i = 1, 2, . . ., n. This rotation represents the "attitude" of the new position of 
the spacecraft in the coordinate system defined by the old position of the space- 
craft control coordinate system c. The following "before" and "after" figures 
may be instructive. 
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From these figures and matrix analysis, it can be seen that the rotation R 
which satisfies R v = w has the form 

/ Cdi p -Sdi/» 

H = 

y Sdt p Cd \p 

This represents an "attitude” with "right ascension" -dp. 



An SCPS computer program - DOAOP - is used by MESS to compute the 
rotation matrix M which minimizes the function 


f (M) = 


L 

i = 1 



for given input vectors v^ , w £ , i = 1, 2, . . n. Thus DOAOP can be used to 
find the rotation matrix which "best" (in the sense of least squares) estimates 
the attitude change required above. For small pitch and yaw misalignments 
66, dip, resp., the attitude matrix returned by DOAOP will have the nominal 
form 

( 1 -dip 66 

dip 10 

-69 0 1 



4.5. Analysis for t he IRU Exercise. The IRU is designed to provide precise 
slew sensing about each of the spacecraft control coordinate axes in either the 
positive or negative direction (taken in the right-hand sense). First-order 
deviations in slew-axis alignment are provided in the SCPS spacecraft model 
by the matrix 
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i + an = 



The nominal values of the off-diagonal elements do;.., i ^ j, in the above matrix 
are all zeroes. The i th column (i = 1, 2, 3) of the above matrix is a first-order 
unit vector giving the actual (misaligned) i th slew axis (1 = roll, 2 = pitch, 3 = 
yaw) represented in the control coordinate system. 

First-order deviations in slew-angle sensing are provided in the SCPS 
spacecraft model by slew angle scale factor calibration parameters. There are 
six such, one each for the positive and negative directions of slew about the 
three slew axes. These parameters are taken to be small corrections £ which 
have to be added to the nominal scale factor 1 to produce the correct scale factor 
for the slew angle. That is, if a slew of angle 6 is commanded about the i-axis 
(1 = roll, 2 = pitch, 3 = yaw) in the j-direction (1 = positive, 2 = negative), the 
slew which actually occurs has angle (1 + £. . ) 9 . 

The approach taken in the IRU Exercise analysis is to construct a system 
of linear equations which relate the small unknown changes in the values of the 
alignment and calibration parameters, to the small errors in terminal slew 
attitude which result from those changes. If a large number of' such equations 
is available, the unknown changes in parameter values may be estimated by 
standard linear estimation techniques. 


Symbolically, each error e measured in a certain terminal slew attitude 
may be written as a linearized function of these unknown parameter changes A^: 
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e= ZZ A# 

i 


The coefficients C^, which vary from observation to observation, may be collected 
into a row vector of coefficients, Cj. The unknown changes A g in parameter 
values, which are assumed to remain the same from observation to observation, 
may be collected into a column vector of parameters, P: 


P = 


Then the error e may be expressed 
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If several errors e are available, from the same or other slews, one can 
write an equation like the above for each. The row of coefficients, £ J, will vary 
as various errors are used, but the column of unknowns, P, is assumed to re- 
main fixed. Now a large number of error observations from many different 
slews may be collected into a column array, E, and the rows associated with 
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these errors can be assembled into a coefficient matrix C. Then we can write 
the following, called the IRU Estimation Equation: 

CP = E 

The above equation will be constructed operationally as follows. After the 
errors e f , e p , e y in roll, pitch and yaw for a certain slew are read, the 36 co- 
efficients C r£ ’ C p<f ’ C y£ which relate these errors to the 12 unknown parameter 
changes A g are calculated. (Formulas for calculating these coefficients are 
derived later in this paper.) These data are then assembled into the next three 
rows of the Estimation Equation, as shown diagrammatically in Figure 7. 


column 



Figure 7. Data Placement in IRU Estimation Equation 

The design philosophy of the IRU Exercise then is to provide a large number 
of slews of a wide range of angles about all three axes, with an error readout 
at the end of each slew. In this way sufficient error terms E can be collected 
and coefficients £ computed for the Estimation Equation, to enable a valid 
statistical solution for the parameters P to be determined. The remainder of 
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this section will be devoted to developing formulas for computing the coefficients 


e. 


Suppose a slew is commanded, say a three-legged slew, from an initial 
stellar-referenced attitude, and the slew is allowed to take place under IRU 
control. Suppose that S x , S 2 , S 3 represent the nominal coordinate transforms 
resulting from the first, second, third legs, respectively. Suppose further that 
Sj , S 2 , S^ represent the actual coordinate transforms resulting from the mis- 
aligned slew legs. Then we have the coordinate system flow diagram shown in 
Figure 8. 


initial 

attitude 



In Figure 8, the transformation S from a 2 to a 2 is given by 

O — C # c* cT oT qT 
^"^ 3^2 ^2 ^3 

The attitude a 2 is by definition the attitude in which the coordinate system which 
is stellar referenced at a 2 actually ended up. This coordinate reference was 
nominally to have ended up at a 2 . Therefore if some combination of stellar- 
referenced sensors (GST, BST, FES) is allowed to control the spacecraft sub- 
sequent to settling out (at attitude a') from the IRU-controlled slew, the space- 
craft will be brought (except for misalignment of the stellar reference) to 
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attitude a 2 . When this happens, the IRU will take on roll, pitch and yaw errors 
e r , e p , e y . These must therefore be the angles which determine the small-angle 
rotation S: 


Q _ c' Q' C* qT oT eT _ 
^ ~ ^2 ^3 “ 




e 

y 


1 


-e 


r 



Each error e , e , e can be considered a function of the six nonorthogonality 
r p y 

parameters do) mn (m, n = 1, 2, 3; m/n) and the six slew calibration parameters 
(p = 1, 2, 3; q = 1, 2). We can write formally e s = e & (d« 12 , do) 13 , . . ., £ 32 ), 
s = r, p, y. Since each e s has nominal value e s = 0, we can relate the errors e s 
to changes A dc^ , A in the parameters. To first order, 
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The above equation relates the changes Ado> 12 , Ada> 13 , ...» A£ 32 in the param- 
eters of interest to errors e r , e p , e y measured in the operating environment. 
The coefficients 



are the desired coefficients C s ^ of the IRU Estimation Equation. 

Let § be any of the 12 parameters to be estimated. Then we have 
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■jj (s; sj s; sj sj sj> 


as: as: as; 

£ Q' O' cT oT cT e # £ C' eT oT cT , o # C # L qT oT oT 

3^: S 2 b l °1 & 2 b 3 + & 3 3^r b l & 1 & 2 S 3 + 5 3 S 2 3^: *2 ^3 


Evaluating this equation at the nominal state 0, S. = S!, i = 1, 2, 3, yields the 
desired coefficients 



r, p, y 


which appear as elements of the array 
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Hence the coefficients C s ^ depend on the slew S and the partials 



= 1, 2, 3 


The subscripts i,0 will be dropped hereinafter for readability. 

H £=<•"«, for some m, n, we will have B S'/d'g = 0 unless the slew S’ takes 
place about the n-axis. If the slew S' does take place about the n-axis, it is 
shown in the Appendix that 

J^L- = (1 -C9) (e m e T n + S n ej) - S 9 e m ; m, n = 1, 2, 3; m * n 

mn 

where 9 is the angle slewed, including algebraic sign; e m , e n are the Kronecker 

r\j 

basis vectors along the m, n axes, respectively; e m is the slew-symmetric ma- 
trix associated with the vector e m according to the correspondence 



If § - £ pq for some p, q, we will have B S’/B£ = 0 unless the slew S f takes 
place about the p-axis in the q-direction. In this latter case, suppose a slew S 
of angle 9 is commanded about the p-axis in the q-direction. Then the slew S’ 
which will actually occur is (1 + £ pq ) 9, and it is shown in the Appendix that 

= -6 (c 9 e p + S9(I - ; p aj>) ; p = 1, 2, 3; q - 1, 2 

^pq 
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6. Appendix: Derivation of Formulas for Computing 

It is well known that the matrix S of a rotation about an axis represented by 
the unit vector 



through the angle 9 can be expressed in "axis-angle" form: 


s = ce 



where e is the skew- symmetric matrix associated with the vector e according 
to the correspondence 


— x. 


e = 


-x. 


r 

X 3 

V- 

To compute d S'/3£ when £ = d first observe that each misaligned co- 
ordinate vector e’ is represented in the nominal coordinate system by the n th 
column of the array I + dC, where dfl is the array of nonorthogonality parameters: 
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Then a rotation about the misaligned axis e n through angle 6 has the matrix 
representation 


S" = C <9 I + (1 = C d) e* e* T - S 0 e' 


= ce i + (i -c£) (i + dfi) n (i + d D)i - se (T + u~i) n 

= Cdl + (l-C< 9 )(I n lJ + dn n lT + l n dfiJ)-Si 9(I n+ d\) 

including terms through first order. Now only the parameters ^ 0J mn , m 4 n, 
appear in dft n . For these parameters, 

3 d a 

n 

3dw ” Sm 


and also we notice that I n is just e n . Hence 


3 S' 


d d co 


mn 


— - — (VdQ i T + 1 dn T > (i - ce) - se dh A 

ddco V n n n V ' "/ 

mn 



+ e e T ) (1 - C 9) - S 9 e 

n m' ' ' m 


To compute 3S’/3^ when^ = £ , again the axis-angle form of a slew is 

utilized. Now a slew S’ which takes place about an axis represented by the unit 
vector e p in the q-direction through an angle (1 + £ pq ) 6 has the form 

S' =C(1 + i vo )91 .+ (1 - cos (l + £ pq )S) e p e£ - S (1 + y^e p 
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Differentiating yields 

3S # 


-0S(1+'C )0I+0S(l+£ )9e e T - 9 C (1 + £ )9e 
v pq' v ^pq y p p v b pq y p 


pq 


,--d(c(l + i ) 9e + S(I + £ ) 9 (I - e e T )^ 

\ v 5 pq' P v ^pq y v p P ') 

Since the derivative is needed only to zero order, this may be evaluated when 


£ = 0 : 
5 pq 


r 3 S' 

i s pq 


A = - 6 (c9 e p + S0 (I - e p ej)} . 

/ O 
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